Effect of charged particle's multiplicity fluctuations on flow harmonics in 

even-by-event hydrodynamics 
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In nucleon-nucleon collisions, charged particle's multiplicity fluctuates. We have studied the 
effect of multiplicity fluctuation on flow harmonics in nucleus-nucleus collision in event-by-event 
hydrodynamics. Assuming that the charged particle's multiplicity fluctuations are governed by the 
negative binomial distribution, the Monte-Carlo Glauber model of initial condition is generalised 
to include the fluctuations. Explicit simulations with the generalised Monte-Carlo Glauber model 
initial conditions indicate that the multiplicity fluctuations do not have large effect on the flow 
harmonics. 

PACS numbers: 47.75.+I, 25.75.-q, 25.75.Ld 



I. INTRODUCTION 

In recent years, there is much interest in event-by- 
event hydrodynamics. Event-by-event hydrodynamics 
takes into account that in nucleus- nucleus collisions, par- 
ticipant nucleon positions fluctuates from event to event. 
Effect of participant nucleon position fluctuations is most 
prominent on the azimuthal distribution of the produced 
particles. In a non-zero impact parameter collision be- 
tween two identical nuclei, the collision zone is asym- 
metric. Multiple collisions transform the initial asymme- 
try into momentum anisotropy. Momentum anisotropy 
is best studied by decomposing it in a Fourier series, 
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(1) 

cj> is the azimuthal angle of the detected particle and tp n 
is the plane of the symmetry of initial collision zone. For 
smooth initial matter distribution (obtained from geo- 
metric overlap of density distributions of colliding nu- 
clei), plane of symmetry of the collision zone coincides 
with the reaction plane (the plane containing the impact 
parameter and the beam axis), tp n = ^up,^n. The odd 
Fourier coefficients are zero by symmetry. However, fluc- 
tuations in the positions of the participating nuclcons 
can lead to non-smooth density distribution, which will 
fluctuate on event- by- event basis. The participating nu- 
cleons then determine the symmetry plane (V'pp), which 
fluctuate around the reaction plane [l|. As a result odd 
harmonics, which were exactly zero for smoothed initial 
distribution, can be developed. It has been conjectured 
that third harmonic U3, which is response of the initial 
triangularity of the medium, is responsible for the ob- 
served structures in two particle correlation in Au+Au 



collisions @,0,0,[l],[i],0. The ridge structure in pp 
collisions also has a natural explanation if odd harmonic 
flows develop. Recently, ALICE collaboration has ob- 
served odd harmonic flows in Pb+Pb collisions [8]. In 
most central collisions, the elliptic flow (i^) and trian- 
gular flow (^3) are of similar magnitude. In peripheral 
collisions however, elliptic flow dominates. 

It is well known from experiments that in nucleon- 
nucleon collisions, charged particles multiplicity fluctu- 
ates strongly. Experimental data in limited phase space 
as well as in full phase space are well described by the neg- 
ative binomial distribution Recently, in [ll|, in 
a model calculation with Color Glass Condensate (CGC) 
initial energy density distribution, it was observed that 
the experimental multiplicity distribution in d+Au colli- 
sions at RHIC are better explained if multiplicity fluctua- 
tions in NN collisions are included. It was also shown that 
various moments of the eccentricity of the collision zone 
in nucleus- nucleus collisions are affected by the multiplic- 
ity fluctuations in NN collisions. Fluctuation dominated 
moments, e.g. triangularity, increase substantially when 
the effect of multiplicity fluctuations is included. Effect 
of multiplicity fluctuations on flow coefficients however 
is not studied in [TTJ . If, higher moments of the eccen- 
tricity of the collision zone increase substantially due to 
multiplicity fluctuations, an increase in flow harmonics is 
expected. 

In general, event-by-event hydrodynamics do not ac- 
count for the multiplicity fluctuations in NN collisions. 
In the present paper, in a simple model, we incorporate 
the effect of multiplicity fluctuations in event-by-event 
hydrodynamics and study its effect on the flow harmon- 
ics. The model is substantially different from that used 
in While in [TO], Color Glass Condensate (CGC) 
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model is used to obtain the initial energy density distri- 
bution, in the present model, we have used Monte-Carlo 
Glauber model to obtain the initial energy density dis- 
tribution, event-by-event. Explicit simulations indicate 
that multiplicity fluctuations do not have large effect on 
the flow harmonics. 
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II. MONTE-CARLO GLAUBER MODEL WITH 
MULTIPLICITY FLUCTUATIONS 



Npart, and (ii) a hard component, which is proportional 
to the number of binary collisions N co u , 



In theoretical simulations of event-by-event hydrody- 
namics, one generally uses the Monte-Carlo Glauber 
model to obtain the initial energy density distribution 
in an event. Details of the Monte-Carlo Glauber model 
can be found in [ijj]. In a Monte-Carlo Glauber model, 
according to the density distribution of the colliding nu- 
clei, two nucleons are randomly chosen. If the trans- 
verse separation between them is below a certain dis- 
tance, they are assumed to interact. Transverse po- 
sition of the participating nucleons is then known in 
each event. The positions will fluctuate from event-to- 
event. If a particular event has N part participants, partic- 
ipants positions in the transverse plane can be labeled as, 
(xi,yi), (x2,y2)~~(x Npart ,y Npart ). Energy density distri- 
bution in the particular event can be obtained by as- 
suming that each participant deposit energy £o in the 
transverse plane, 



e(x,y) ps £ 



N part 



S(x- 



,y-yi) 



(2) 



Fluid dynamical models require continuous density dis- 
tribution and discrete distribution as in EqJ5] cannot be 
evolved in a hydrodynamical model. To use in a hy- 
drodynamic model, the discrete density distribution has 
to be converted into a smooth energy-density distribu- 
tion. This can be done by smearing the discrete partic- 
ipating nucleon positions by some smoothing function, 
6(x -Xi,y- y t ) g(x - x u y - y it £i, &■■), d being pa- 
rameters of the smoothing function g. 



e(x,y) = e ^ g(x - y,Xi,y - y i: Ci, C2-O (3) 



One generally uses a Gaussian smoothing function. 
However, there can be other choices, e.g. in |l3|, a 
Woods-Saxon distribution function was used for the 
smoothing. In the present simulations, we have used a 
Gaussian distribution 



g gauss (x- Xi,y - yi,a) oc e' 



(4) 



of width a— 0.5 fm. 

As it was mentioned earlier, in experiments, charged 
particle's multiplicity fluctuates in NN collisions. How we 
incorporate the effect of multiplicity fluctuations within 
the framework of Monte-Carlo Glauber model? Gen- 
esis of Eq[3] gives a clue. It was observed that in 
nucleus-nucleus collisions, charged particles rapidity den- 
sity {dN c h/dr]) have two components, (i) a soft compo- 
nent which is proportional to the number of participants 



drj 



(1-x) 



Npart dN% 

2 dri 



xN, 



coll " 



dr] 



(5) 



The hard scattering fraction x is small (~ 0.1-0.2) in 
RHIC and LHC energy collisions. If we assume that EqE] 
is valid in all the transverse positions r±, and further 
assume that initial transverse energy density is propor- 
tional to charged particles multiplicity, one can write the 
transverse energy density distribution, in nucleus-nucleus 
as, 



e{x,y) = s 



Np ar t (X,y) 



+ xN co u(x,y) 



(6) 



A large number of experimental data in y/s NN —200 
GeV Au+Au collisions are well explained in (smooth) 
hydrodynamic model with initial energy density as given 
in EqlHl One understand that the Monte-Carlo Glauber 
model initial energy density distribution in Eq|3]is essen- 
tially a generalisation of Eq|6] with the hard scattering 
fraction x=0 and £0 in Eqj3] is essentially a measure of 
charged particles multiplicity in NN collisions. 

Charged particle's multiplicity fluctuations in NN col- 
lisions imply that even in a single event, all the N par t par- 
ticipant pairs do not produce the same multiplicity i.e., 
£0 fluctuates. Considering that in NN collisions, multi- 
plicity fluctuations follow the negative binomial distribu- 
tion, Monte-Carlo Glauber model of initial condition can 
be generalised to include multiplicity fluctuations as, 



e(x,y) = Af 



part 



1=1 



](x~y,x l ,y-y i Xi,C2---), (7) 



In Eqlll (xi,yt) are the transverse positions of the par- 
ticipant nucleons, g is a smoothing functions with pa- 
rameters Ci- n i is a random number following the two 
parameter negative binomial distribution P(n : (n), k). 



(n),k) 



T(n + k) 
T(n+ l)T(k) 



The parameters (n) and k corresponds to the average of 
the distribution and the width of the distribution. Af in 
Eql7] is an overall normalising factor. It can be fixed to 
reproduce experimental observable, e.g. multiplicity or 
transverse momentum distribution. 

Negative binomial distribution parameters (n) and k 
in EqJS] corresponds to pp collisions multiplicity distri- 
butions. However, it is possible that in nucleus-nucleus 
collisions, pp collision characteristic is changed. For ex- 
ample, in |11| . it was shown that experimental data on 
charged particles multiplicity distribution in d+Au colli- 
sions are better fitted with, 



3 



(a) No Nch flue 



(b) Nch fluc.-l 



k = k pp - min(T A (r±),T B (r±))a 



(9) 



In the following, we have simulated 20-30% Au+Au 
collisions at RHIC energy, ^^^=200 GeV. In 
1 /s Arjv =200 GeV Au+Au collisions, average charged par- 
ticles multiplicity per participant — - d ^ h slowly in- 
crease with collision centrality. In 20-30% collision, 
5N art ^t ~ 33 E3- We then nx (n)=3. Simulation 
results for a higher average multiplicity (n}=10 will also 
be shown. For the width parameter, we have considered 
two possibilities, which we call model I and II for N c h 
fluctuations. 



Model - I :k 
Model - II :k 



min(T A (r ± ),T B (r±))a Q (10) 
(11) 



In y / s=540 GeV pp collisions, in central rapidity, 
kpp rsi 2 [loj . kpp decreases logarithmically with energy. 
We then fix k pp =2. The difference between model I and 
II is the following: in model II, in EqJT] is a random 
number following a negative binomial distribution with 
parameters, (n) = 3 and k = k pp = 2. In model I, rij 
is still a random number following negative binomial dis- 
tribution with average (n)=3, but now k is not a fixed 
parameter, rather it depend on the participant nucleon 
positions. Here, I briefly mention the essential difference 
of the present model with that in [llj . In 11], initial en- 



ergy density was obtained from a color glass condensate 
model. In CGC model initial conditions, initial energy 
density is assumed to be proportional to gluon multiplic- 
ity, 



dN A+B^g 

dyd 2 TT 



K-^- j f d*k T a x (Q) (12) 
N G -I J Pt Jt 



|p T + k T | \ /|p T -k T | 

o > x l * n 



where ^(kx,x) is the unintegrated gluon distribution 
function. Above equation is interpreted as average local 
multiplicity. In each cell V 2 tt, of the transverse plane, 
the actual multiplicity is a Negative Binomial distribu- 
tion random variable. In other words, in |Tlj |. m each 
fluid cell of the transverse plane is associated with fluctu- 
ations. In the present model, fluctuations are associated 
with the participant nucleon only. 



III. HYDRODYNAMIC EQUATIONS, 
EQUATION OF STATE, INITIAL CONDITIONS 

With generalised Monte-Carlo Glauber model initial 
condition for initial energy density, space-time evolution 
of the fluid, in each event is obtained by solving the 
energy-momentum conservation equations, 




FIG. 1: (color online) Energy density distribution in a typical 
Monte-Carlo Glauber model event in 20-30% Au+Au colli- 
sions, (a) without any multiplicity fluctuations, (b)-(c) with 
multiplicity fluctuations following model-I and II respectively. 



{e+p)u^u v - Pg h 



d^ u = 0, 



(13) 
(14) 



where e and p are the energy density and pressure re- 
spectively, u is the hydrodynamic 4-velocity. We have 
disregarded any dissipative effect. Assuming boost- 
invariance, hydrodynamic equations are solved in (r — 



= ~ In jzr-) coordinate system, with the 



code AZHYDRO-KOLKATA [14|. Hydrodynamics equa- 
tions are closed with an equation of state (EoS) p = 
p(e). Currently, there is consensus that the confinement- 
deconfinement transition is a cross over. The cross over 
or the pseudo critical temperature for the transition is 
T c sal 70 MeV [l5Wl8|. In the present study, we use an 
equation of state where the Wuppertal-Budapest [l5|, [l?} 
lattice simulations for the deconfined phase is smoothly 
joined at T = T c = 17 '4 MeV, with hadronic reso- 
nance gas EoS comprising all the resonances below mass 
m res ~2.5 GeV. Details of the EoS can be found in [l9j . 

In addition to the initial energy density for which we 
use the Monte-Carlo Glauber model, with or without 
multiplicity fluctuations, solution of hydrodynamic equa- 
tions requires to specify the thermalisation or the initial 
time r, and fluid velocity (v x (x, y), v y (x, y)) at the initial 
time. A freeze-out prescription is also needed to convert 
the information about fluid energy density and velocity 
to invariant particle distribution. We assume that the 
fluid is thermalized at t,= 0.6 fm and the initial fluid ve- 
locity is zero, v x (x,y) = v y (x,y) = 0. The freeze-out 
is fixed at Tp=130 MeV. We use Cooper-Frye formal- 
ism to obtain the invariant particle distribution of 7r~ 
from the freeze-out surface. Resonance production is in- 
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FIG. 2: (color online) Charged particles transverse momen- 
tum distribution with and without multiplicity fluctuations. 



eluded. Considering that pions constitute ~ 20% of all 
the charged particles, 7r _ invariant distribution is mul- 
tipled by the factor 2 x 1.2 to approximate the charged 
particle's invariant distribution. From the invariant dis- 
tribution, harmonic flow coefficients are obtained as pol ] , 



dN 



(15) 



dyp T dpT 



dypT dpT dep 



dN 
dy 



(16) 



In a boost-invariant version of hydrodynamics, flow co- 
efficients are rapidity independent. Present simulations 
are suitable only for central rapidity, y «0, where boost- 
invariance is most justified. Hereafter, we drop the ra- 
pidity dependence. ip n in Eqs ll5U"6l is the participant 
plane angle for the n-th flow harmonic. We characterise 
the asymmetry of the initial collision zone in terms of 
various moments of the eccentricity [Ej] , [H , Q , 



inip n 



J J e{x,y)r n e l2 ' t 'dxdy 
J J s(x,y)r n dxdy 



= 2,3,4,5(17) 



62 and £3 are called eccentricity and triangularity. In 
the following, 64 and £5 will be called rectangularity and 
penta-angularity respectively. They measure the square- 
ness and five sided-ness of the initial distribution respec- 
tively. EqJT7]also determine the participant plane angle 
ipn- 



IV. RESULTS 

To study the effect of multiplicity fluctuations in 
event- by- event hydrodynamics, we have simulated 20- 
30% Au+Au collisions for N event =100 events. Recently 
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FIG. 3: (color online) Simulation results for event averaged 
asymmetry parameters e n , n=2-4. e„ marginally increase 
when multiplicity fluctuations are included. 



in [22J, it was shown that with fluctuating initial condi- 
tions, event averaged as well as variance of elliptic flow 
and triangular flow remain approximately unaltered for 
-^ r el , e „ t =50-2500. N event — 100 is then sufficiently large 
to study the effect of multiplicity fluctuations in event- 
by-event hydrodynamics. In each event, Monte- Carlo 
Glauber model participant positions are smoothed with a 
Gaussian function of width cr=0.5 fm. Multiplicity fluc- 
tuations are accounted either with model I or model II. 
We also simulate events without any multiplicity fluctua- 
tions. The overall normalising factor Af is fixed such that 
event averaged charged particles multiplicity reproduces 
the experimental value, ~ 276 in 20-30% Au+Au 

collisions. 



In FigHJ effect of multiplicity fluctuations on the ini- 
tial energy density distribution, in a typical Monte-Carlo 
event is shown. Panel (a) shows the energy density dis- 
tribution without any multiplicity fluctuations. Energy 
density distributions with multiplicity fluctuations are 
shown in panel (b) and (c) respectively for model-I and 
II. Effect of multiplicity fluctuations is to smooth the en- 
ergy density distribution. The fine structures which are 
seen in panel (a) get diffused when multiplicity fluctua- 
tions are included. 



A. Transverse momentum spectra 

In Fig[2J simulation results for the charged particle's 
transverse momentum distribution are shown. The lines 
represent the event averaged values. The standard errors 
are small and are not shown here. The black lines in 
FigO are the simulated spectra when multiplicity fluc- 
tuations are neglected. The red and green lines are simu- 
lated spectra with multiplicity fluctuations included fol- 
lowing model-I and II respectively. Transverse momen- 
tum distribution is marginally affected due to multiplic- 
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FIG. 4: (color online) Probability distribution of initial spa- 
tial anisotropy P{(-2), P{^s), P( e s)y Pfa) an d ^ > ( e s) i n simu- 
lated events are shown. In each panel, the black colors repre- 
sent simulations without multiplicity fluctuations. Multiplic- 
ity fluctuations with model I and II are indicated by red and 
green colors. 



ity fluctuations. Indeed, in the figure, one can not dis- 
tinguish between spectra obtained without multiplicity 
fluctuations and with multiplicity fluctuations I. For mul- 
tiplicity fluctuations II, the spectrum is marginally hard- 
ened. 



B. Integrated flow coefficients 

In a hydro dynamical model, development of the flow 
coefficients depend on the initial spatial anisotropy pa- 
rameters (e n ) characterising the collision zone. Approx- 
imately, larger the initial asymmetry, more is the flow. 
Does the initial asymmetry parameters e„ changes when 
multiplicity fluctuations are included in the Monte- Carlo 
Glauber model of initial condition? Simulation results 
for the event averaged asymmetry parameters e n and 
its standard errors are shown in Fig|3] (e„) increase 
marginally when multiplicity fluctuations are included 
for. The increase is barely more in model-II than in 
model I. Marginal effect of multiplicity fluctuations on 
asymmetry parameter is also evident in Fig 01 where, in 
four panels, histograms for the probability distribution 
of the initial eccentricity (£2), triangularity (£3), rectan- 
gularly (£4) and penta-angularity (£5), in the simulated 
events are shown. The black bars correspond to simu- 
lations without multiplicity fluctuations. The red and 
green bars are obtained when multiplicity fluctuations 
are included according to model I and II respectively. 
Qualitatively, probability distribution of asymmetry pa- 
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FIG. 5: Simulated flow coefficients v n , n=2,3,4 and 5 with 
and without energy fluctuations. 



rameter e n remain similar. The result is not unexpected. 
In Monte-Carlo Glauber model, the asymmetry parame- 
ters are determined mainly by the positions of the partic- 
ipant nucleons (see Eq |17[) , weighted by the energy den- 
sity. Since the participant positions remain unaltered the 
change in e n due to multiplicity fluctuations is small. 

In FigO simulated (integrated) flows, v n , n—2-5, with 
and without the multiplicity fluctuations are shown. 
The symbols represent the event averaged values, the 
bars the standard errors. Event averaged flow coeffi- 
cients marginally increase due to multiplicity fluctua- 
tions. However, the increase is small. In model I of 
fluctuations, flow coefficients v n , n=2-5 increase by less 
than 5%. In model II also, the increase is modest, less 
than 5% in v 2 and V3. In higher flow coefficients V4 and 
V5, the increase is marginally more, ~10-15%. The result 
is in agreement with our observations that initial asym- 
metry parameters depend marginally on the multiplicity 
fluctuations in NN collisions. 

Recently, in [23| . correlation between elliptic (1*2) and 
triangular (^3) flow with initial eccentricity (£2) and ini- 
tial triangularity (£3) was studied. It was shown that in 
event- by- event hydrodynamics, elliptic flow is strongly 
correlated with initial eccentricity. Comparatively weak 
correlation was observed between triangular flow and ini- 
tial triangularity. Correlation between higher flow coef- 
ficients, U4, V5 and initial spatial asymmetry parameter 
£4, £5 were not studied. It is interesting to study the 
correlation between higher flow coefficients V4, v$ and 
asymmetry parameters £4 and £5. Effect of multiplicity 
fluctuations on the correlation between different flow co- 
efficients and asymmetry measures, is also interesting. In 
Fig|6l simulated flow coefficients v n are plotted against 
the initial asymmetry parameters £„. Simulation results 
with and without multiplicity fluctuations are shown in 
separate panels. One observes that elliptic flow v-i is 
strongly correlated with initial eccentricity £2. The tri- 
angular flow V3 also appear to be correlated with initial 
triangularity, though the degree of correlation is less than 
in elliptic flow. However, correlation between higher flow 
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FIG. 7: (color online) Simulated V2, vz, i>4 and «5 are shown 
in four panels. In each panel, flow coefficients, without mul- 
tiplicity fluctuations (the black line), with multiplicity fluc- 
tuations with model I (the red line) and with model II (the 
green line) are shown. Average (n) of the negative binomial 
distribution is fixed at {n)—3. The blue circles are simulation 
results with model II of fluctuations but with (n)=10. 



FIG. 6: Effect of multiplicity fluctuations on the correlation 
between flow coefficients and initial asymmetry parameters, 
(«2,£2) («3,£3), ("4,£4) and (v 5 ,e 5 ), are shown. 

coefficients w 4 and v 5 with spatial asymmetry e 4 and e§ 
appears to be weak. Effect of multiplicity fluctuations 
on the correlation appear to be marginal again. In 23] a 
quantitative measure was defined to quantify the corre- 
lation between flow coefficients and initial spatial asym- 
metry measure. A modified form is used here to measure 
the correlation, 



TABLE I: Correlation measure for flow harmonics v n , n=2-5 
in the event-by-event hydrodynamics with and without mul- 
tiplicity fluctuations. 
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measure 






V2 




V 4 


Vs 


no Nch flue. 


0.980 


0.873 


0.633 


0.712 


Nch flue. -I 


0.974 


0.865 


0.538 


0.727 


Nch fluc.-II 


0.971 


0.869 


0.466 


0.689 



„ ( \ — 1 _ Z-iiV V n\ e n) v n,st.line\(-n)\ (-\s\ 

Cmeasure essentially measures the dispersion of the 
simulated flow coefficients from the best fitted straight 
line, relative to completely random flow coefficients. It 
varies between and 1. If flow coefficients are perfectly 
correlated then v n oc e n and C meaS ure is identically unity. 
For completely random flow coefficients, C meastire =0. To 
obtain an even ground for Comparison Of Cmeasure 

for dif- 
ferent flow coefficients, the flow coefficients (v n ) and the 



asymmetry parameters (e„) are scaled to vary between 
and 1. In tableU we have listed the C measure values 
obtained in the present simulations. 

Let us first discuss the correlation measures without 
any multiplicity fluctuations. Elliptic flow is strongly 
correlated with initial eccentricity, C measur e(v2) ~ 0.98. 
Correlation in triangular flow is comparatively weak, 
C m easure{vd) ~ 0.87. Correlation between higher flow 
harmonics, U4 and the initial asymmetry parameter £4 
or between V5 and £5 is much more weaker than that 
in elliptic or triangular flow, C measU r e (v4,) ~ 0.63 and 
Cmeasure{vh) ~ 0.71. If departure of C meaS ure from unity 
is a measure of flow uncorrelated with the initial eccen- 
tricity measure, only ~ 2% of elliptic flow is uncorrelated 
with initial eccentricity. Uncorrelated flow is more ~13% 
for triangular flow. For higher harmonics, U4 and V5, un- 
correlated flow is much more, ~ 40% in U4 and 30% in V5. 
Multiplicity fluctuations appear to reduce the correlation 
between flow coefficients and asymmetry measures. How- 
ever, with the exception of D4, correlation between flow 
coefficients and asymmetry measures are reduced only 
marginally (less than a few percent). Only in u 4 , corre- 
lation is significantly reduced, e.g. by ~ 15% in model I 
of fluctuations and by ~ 26% in model II of fluctuations. 
We conclude that with the exception of flow harmonic t>4, 
the correlation between flow coefficients and initial asym- 
metry parameter remain largely unaffected by the inclu- 
sion of multiplicity fluctuations in NN collisions. Why 
the correlation between flow harmonic W4 and the asym- 
metry measure €4 is particularly sensitive to multiplicity 
fluctuations is not understood. More study is required to 
resolve the issue. 
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C. Differential flow coefficients 

As it is with the integrated flow coefficients, differen- 
tial flow coefficients also depend marginally on multiplic- 
ity fluctuations. In FigjTl in four panels, simulation re- 
sults for the differential flow coefficients, V2(pt), v^(pt) : 
va{pt) and ws(pt) & re shown. Elliptic flow is hardly 
modified by inclusion of fluctuations. Detailed inspec- 
tion indicate that in the pr range 1-2 GeV, differential 
elliptic flow is decreased by less than when multi- 

plicity fluctuations are included. In triangular flow, flow 
is marginally increased due to multiplicity fluctuations. 
In the pt range 1-2 GeV, in model I of fluctuations, the 
increase is about ~4-5%. The increase is marginally more 
^8% in model II. Higher flow coefficients V4 and ^5 also 
increase marginally, e.g. by 6-7% and by 10-12% respec- 
tively in the pr range 1-2 GeV. We conclude that differ- 
ential flow coefficients V2{pt)-V5(pt) are only marginally 
affected by multiplicity fluctuations in NN collisions. 

Negative binomial distribution has two parameters, the 
width k and average (n). In the simulations presented 
here, two possibilities for the width parameter k was con- 
sidered. The average of the distribution was kept fixed at 
(n)=3. To check whether or not simulated flows depend 
sensitively on the average of the distribution, we have 
simulated 100 events with (n)=10, k = k pp = 2. The 
results (the blue circles) are shown in FigJTJ Simulated 
flows with (n)=10 can hardly be distinguished from the 
simulated flows with (n)=3. 

V. SUMMARY AND CONCLUSIONS 

To summarise, we have studied the effect of multiplic- 
ity fluctuations in NN collisions on flow harmonics in 



^^=200 GeV Au+Au collisions. The Monte-Carlo- 
Glauber model for the initial transverse energy distri- 
bution is generalised to include charged particle's mul- 
tiplicity fluctuations in NN collisions. Multiplicity fluc- 
tuations are assumed to be governed by the negative bi- 
nomial distribution, with two parameters, the average 
(n) and the width k. With the generalised Monte-Carlo- 
Glauber model initial transverse energy distribution, we 
have simulated 20-30% Au+Au collisions event-by-event. 
Two possibilities for the width parameter k are consid- 
ered, (i) model I: k = k pp ■ min(Tji(r±),TB{r±))ijQ, i.e. 
the width depend on the local position of the participant 
nucleons, and (ii) model II:fc — k pp , i.e. the width remain 
same as in pp collision. In the simulations, we use k pp =2. 
The average the distribution is kept fixed at (n)=3. The 
simulations are constrained to reproduce experimentally 
charged particles multiplicity in 20-30% collisions. Ex- 
plicit simulations indicate that in event-by-event hydro- 
dynamics, multiplicity fluctuations do not play signifi- 
cant role in the development of the flow harmonics. In- 
tegrated flows are marginally increased with multiplic- 
ity fluctuations. Increase is barely more in model II of 
fluctuations than in model I. The differential flow is also 
increased marginally. Indeed, within the uncertainties in- 
tegrated and differential flow coefficients remain largely 
unaltered irrespective of multiplicity fluctuations. We 
have also studied the correlation between different flow 
harmonics and initial asymmetry measures. With the 
exception of fourth flow harmonic (v^), correlation be- 
tween flow coefficients and asymmetry measures also re- 
main largely unaltered. We conclude that in Monte-Carlo 
Glauber model of initial condition, flow harmonics are 
not largely affected by the charged particle's multiplicity 
fluctuations in NN collisions. 



[1] S. Manly et al. [PHOBOS Collaboration], Nucl. Phys. A [12 
774, 523 (2006) 

[2] A. P. Mishra, R. K. Mohapatra, P. S. Saumia, A. M. Sri- [13 

vastava, Phys. Rev. C81, 034903 (2010). 
[3] A. P. Mishra, R. K. Mohapatra, P. S. Saumia, A. M. Sri- [14 

vastava, Phys. Rev. C77, 064902 (2008). [15 
[4] J. Takahashi, B. M. Tavares, W. L. Qian, R. Andrade, 

F. Grassi, Y. Hama, T. Kodama, N. Xu, Phys. Rev. Lett. [16 

103, 242301 (2009). 
[5] B. Alver, G. Roland, Phys. Rev. C81, 054905 (2010). [17 
[6] B. H. Alver, C. Gombeaud, M. Luzum, J. -Y. Ollitrault, 

Phys. Rev. C82, 034913 (2010). [18 
[7] D. Teaney, L. Yan, Phys. Rev. C83, 064904 (2011). [19 
[8] [ ALICE Collaboration ], Phys. Rev. Lett. 107, 032301 

(2011). [20 
[9] G. J. Alner et al. [UA5 Collaboration], Phys. Lett. B 160, [21 

199 (1985). 

[10] G. J. Alner et al. [UA5 Collaboration], Phys. Lett. B 160, [22 
193 (1985). [23 

[11] A. D umitru and Y . Nara, Phys. Rev. C 85, 034907 (2012) 
|arXiv:1201.6382l [nucl-th]] . 



B. Alver, M. Baker, C. Loizides and P. Steinberg, 
larXiv:0805.44i"T1 [nucl-ex]. 

M. .Rihan Haque, V. Roy and A. K. Chaudhuri, Phys. 

Rev. C 86, 037901 (2012) 

A. K. Chaudhuri. larXiv:0801.3"T80l [nucl-th] . 

Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz and K. K. Sz- 

abo, Nature 443, 675 (2006) 

Y. Aoki, S. Borsanyi, S. Durr, Z. Fodor, S. D. Katz, 
S. Krieg and K. K. Szabo, JHEP 0906, 088 (2009) 
S. Borsanyi e t al, JHEP 1011, 077 (2010) 
[arXiv:1007.2580l [hep-lat]] . 

Z. Fodor, J. Phys. Conf. Ser. 230 (2010) 012013. 

V. Roy and A. K. Chaudhuri, Phys. Lett. B 703, 313 

(2011) 

Z. Qiu and U. W. Heinz, Phys. Rev. C 84, 024911 (2011) 
S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 
71, 034908 (2005) [Erratum-ibid. C 71, 049901 (2005)] 
A. K. Chaudhuri, Phys. Lett. B 710, 339 (2012) 
A. K. Chaudhuri, Phys. Lett. B 713, 91 (2012) 



